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By W. Thomas Blackshear and Donald G. Eide 
Langley Research Center 


SUMMARY 


The equilibrium free surface of a liquid contained in a tank which is 
being simultaneously subjected to a translational acceleration and a rotation 
about an axis fixed with respect to the tank has been investigated. A general 
differential equation has been derived which governs the three-dimensional equi- 
librium free surface as a function of an inertial field due to translational 
acceleration, angular velocity, and surface tension. Because of the nonlinear 
character of this equation, a complete closed form solution is not presently 
possible. However, a closed form solution is presented which is valid when the 
capillary forces acting on the liquid may be neglected. 

When the capillary forces are of the same order of magnitude as the iner- 
tial forces acting on the liquid, attention has been restricted to two dimensions 
and solutions were obtained by numerical integration in an effort to display the 
effects of various parameters on the free surface shape. Representative two- 
dimensional solution curves are presented and discussed for various magnitudes 
of translational acceleration, rotational rate, liquid-to-solid contact angle, 
and vapor volume. By obtaining a lattice of two-dimensional surfaces, the pos- 
sibility of determining a three-dimensional surface is indicated. 


INTRODUCTION 


Liquid fuel systems and fluid subsystems for space vehicles must be 
designed for an environment where the forces acting on the liquid due to trans- 
lational and rotational accelerations are of the same order of magnitude as the 
capillary forces acting on the liquid. When a liquid free surface exists, the 
motions and eventual equilibrium position of the liquid will be' determined by 
these small forces during some part of almost any space mission. 

A knowledge of the liquid equilibrium configuration in a low force environ- 
ment will be needed in order to design systems such as those requiring restart 
of liquid fuel rockets after periods of coasting flight, evaporative temperature 
control devices, and effective liquid fuel tank venting. The determination of 
liquid equilibrium configurations is also a necessary first step in the applica- 
tion of pertubation analysis to the problem of fuel sloshing. 



Although some experimental and analytical work hearing on some aspects of 
the problem (refs. 1 to 6) has been previously carried out, none of these 
investigations have included the effect of rotations about an axis perpendicular 
to an axis of symmetry of the liquid and removed from the surface of the liquid. 

The purpose of this investigation is to determine the equilibrium free sur- 
face shape of a contained liquid experiencing inertial forces due to thrust and 
rotation that are of the same order of magnitude as the capillary forces acting 
on the liquid. A differential equation is derived which governs the three- 
dimensional free surface of the liquid as a function of an inertial field, 
angular velocity, and surface tension. Under certain restricted conditions 
this equation can be solved and the surface determined in closed form. These 
cases are discussed. In an effort to display the effects of various parameters 
on the free surface shape, attention is then restricted to two dimensions. The 
resulting equation is solved by numerical integration subject to constraints of 
an invariant contact angle at the interface of the liquid, liquid vapor, and 
tank wall, and an invariant vapor area. Representative solution curves are 
shown for various force loadings. In the solution of the equation an unknown 
constant is required In order that the solution satisfy all boundary conditions. 
The constant must be obtained by iterative trial -and-error techniques on an 
electronic computer. An empirical equation for this constant is presented for 
the range of interest of the equations. A method for obtaining a three- 
dimensional free surface by the superposition of a lattice of normal two- 
dimensional surfaces is indicated. 


SYMBOLS 


A 

a 


a x> a y* a z 


_d_ 

dt 

F 


constant defined in equation (ll) 
acceleration due to thrust 

components of acceleration relative to tank-fixed coordinates 
total time derivative 

body forces acting on fluid particle 


-r» 


K 

N Bo 

Ncu 


base vectors for x,y, z coordinate system 
mean curvature of liquid free surface 


Bond number, 


centrifugal number. 


PC0 2 r3 

2a 


2 



p 

Pv 

q. 

R 

r 


V 

X,Y,Z 


press lire in fluid 
pressure in vapor 
velocity of fluid particle 
half -width of tank 
position vector 
vapor area ratio 
inertial coordinates 
tank -fixed coordinates 
liquid-to-solid contact angle 


P 

a 

a 

(JO 

V 


nondimens ionali zed tank-fixed coordinates 
fluid density 

surface-tension coefficient of fluid 
scalar potential of body force 
angular velocity 

components of angular velocity relative to inertial axes 

vector differential operator, ^i+JLj+JLk 

ox ay Sz 

differentiation with respect to inertially fixed coordinates 
differentiation with respect to tank-fixed coordinates 


( ) X ( ) vector cross product 
— > — > 

( ) • ( ) vector dot product 


Subscripts : 

o initial value 


n 


nth value 


Primes denote partial differentiation with £ . An asterisk denotes coordi- 
nate transformation. Arrows over symbols indicate vectors. Dots over symbols 
indicate derivatives with respect to time. 


DERIVATION OF FREE SURFACE EQUATION 


A tank partially filled with an incompressible fluid is subjected to a 
thrust which is constant in magnitude and the direction of which is constant 
with respect to the tank. This tank is 
also experiencing a constant rate of 
rotation about some axis the position 
of which with respect to the tank 
remains fixed. Two systems of coordi- 
nates centered on the axis of rotation 
are defined, one inertial and one fixed 
with respect to the tank. These coor- 
dinate systems are shown in figure 1. 

The equation of motion of a par- 
ticle of the fluid may be written as 
(ref. 7): 



a?| _ ^ vp 

atlj. ‘ F ‘ T 


(i) 


where 

\j denotes differentiation with 
respect to inertially fixed 
coordinates 

q* particle velocity relative to 

inertially fixed coordinates 

— > 

F body forces 

p fluid pressure 

p fluid density 


(a) Three-dimensional view. 



(b) Two-dimensional view. 

Figure 1.- Liquid fuel tank under the 
influence of a constant thrust in 
the negative z -direction and a con- 
stant rate of rotation about the 
Y-y axis. 


The acceleration of the particle in the tank-fixed coordinate system is 
(from ref. 8): 


dq 

dt 


T 



— > — s 

CD X q 


( 2 ) 


where denotes differentiation with respect to tank-fixed coordinates. 

T 


b 



In order to find the equilibrium surface, the rate of rotation of the tank 
is considered to he constant and the motion of the fluid relative to the tank 
is zero. Therefore, equation (2) is simply 


dq* 

dt 


I 


CO X 


(<j? x r ) 



( — > —A 

\cv X r ) • 

la) X r/ 


2 


(3) 


While the tank rotates, it is subjected to a thrust which is constant with 
respect to the tank. The thrust is assumed to create a conservative force 
field a* which may he represented as 

F - -Vft 

so that 

fi = -a* • ? (4) 

Substituting equations ( 3 ) and (k) into equation (l) yields: 


a* • rU 

v a? x r) • 

(a?x 

A 

_ P 

a • j- ~r 

2 



P_ 


- 0 


where 0 is the null vector. 

Integration of equation ( 5 ) yields 


a 


r + 


\a)Xrj 


= c + ? 


( 5 ) 


(6) 


where C is a constant of integration. 

At the liquid-vapor interface, a pressure discontinuity exists. This 
pressure difference is represented by the Young-LaPlace equation (ref. 2) as 


P - P v + aK 


The substitution of this value of p into equation (6) yields the following 
equation for the equilibrium free surface: 


a • 


r + 


cd x r 1 



h: 

p 


The mean curvature of the free surface is therefore 



r + 


|5?x 





( 7 ) 
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Equation ( 7 ) is the most general equation for the equilibrium free surface 
of a liquid in a partially filled container subjected to a combined rotation 
and a carried acceleration field. It is this equation for which a solution is 
desired. 

The mean curvature of a three-dimensional surface is a nonlinear function 
of partial derivatives of the surface which is of such complexity that, in 
general, a complete closed form solution for the free surface is not possible 
at the present time. In addition, the actual boundary conditions are not com- 
pletely understood. In an attempt to provide some answers as to the effects of 
the combinations of capillary, acceleration, and rotation forces, two cases 
which permit the simplification of equation (6) are investigated. 


SOLUTIONS FOR FREE SURFACE 


Inertial Forces Dominant 

Under the assumption that the inertial force due to thrust and/or the cen- 
trifugal force is sufficiently large to merit disregarding the effects of capil- 
lary forces (such as surface tension), multiplying both sides of equation ( 7 ) 
by a, and setting a equal to zero, the following is obtained: 

(^y 2 + co z 2 )x 2 (co * 2 + co z 2 )y 2 (a ^ 2 + a^ 2 )z 2 

+ + — ai a) xy - a) a) zx 

2 2 2 xy-'xz 

- a) z a>yyz + a^x + a^y + a z z - D =■ 0 (8) 


where 



C 


No generality is lost by selecting the coordinate system such that the y-axis 
lies along the axis of rotation. In this case equation (8) reduces to 


2 2 2 2 
ay x co z 

— 2 — + — 2 — + v + ^y + a z z “ D = 0 


(9) 


By expressing equation (9) in canonical form by the coordinate transforma- 


tions : 


a„ 


x* - x + 


2 oy 

y* * y - sr- - 


Sc 


®y 4ayCO y 2 ^a-yCDy 2 
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(11) 


For the case 


a 


y 


=• 0, equation ( 9 ) reduces to 


X 2 + 


fx 

“V 2 




x + 


V 



which is the equation of a right circular cylinder with radius 

1/2 


^y 2 ^y 4 


and the axis of which is the line xz where x = - and z = - If 

2cry 2oiy. 

the equation of the tank is known in the appropriate coordinate system, either 
equation (lO) or (ll) depending upon the existence of a^ may he solved to 
yield the liquid free surface. An example of the surface is shown in figure 2. 
This solution, however, is subject to the constraint that the vapor volume 
which is the portion of the tank not occupied by the liquid must remain constant. 


Solutions With Significant Capillary Forces 


When the capillary forces are significant, the problem is simplified by 
restricting the investigation to two dimensions. It is further assumed that a? 
is parallel to the z-axis and a? lies along the y-axis so that equation ( 7 ) 
reduces to 




x 2 + 




z - 


p + 

V 


PC 


( 12 ) 


At this point the space coordinates are nondimensionalized by the trans- 
formation: 


x = |R 


z = t}R 


By using the expression for curvature in two dimensions, equation (12) 
becomes 




3/2 


o ^5 £ 

f® R / 2 p) P aR 

* + 1 1 * — 


1 - f p v + p c 


where 


V 


d| 
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The following substitutions are now made: 


n Bo = 


paR 2 

a 


A a - f ( P V + pc ) 


N, 


cu 


p03 2 R 5 

2a 


Where N-g 0 is the familiar Bond number and N cu is hereby defined as centrif- 
ugal number. Therefore 


^75 = Ml 2 + l 2 ) + + A (13) 

[i*(v) a ] ' 

Equation ( 15 ) is second order so that its general solution will contain 
two arbitrary constants. (See ref. 9*) In addition, A contains a constant 
of integration so that altogether three independent boundary conditions are 
required to determine a unique solution for the free surface. Three such con- 
ditions are as follows: 


(l) Because of the symmetry of the forces acting on the free surface about 
the T]-axis, the slope of the surface at £ - 0 is zero; that is, — 

| -0 


- 0 . 


(2) At the interface of three different mediums such as liquid, gas, 
and solid, a contact angle is formed. Experimental work on these contact angles 
indicates a characteristic angle is formed by a given liquid in contact with a 
given tank wall material and surface condition. In reference 2 it is shown that 
this contact angle is invariant with inertial loading. In references 4, 5* 
and 6 the influence of the contact angle on the equilibrium liquid surface is 
shown. Therefore, for this problem various contact angles are assumed by 
requiring that the slope of the free surface at the tank wall be equal to the 

, dri 

cotangent of the contact angle; that is, — - 

d| 6=1 


- -cot 0 


C • 


( 3 ) The vapor area bounded by the tank base, walls, and liquid free sur- 
face must be maintained at a constant value. 


In all cases it was assumed that the vapor area was such that a continuous 
liquid free surface intersecting opposite tank walls would be formed. 

In reference 6, the two-dimensional equilibrium free surface was investi- 
gated in the absence of angular velocity. The free surface equation derived in 
reference 6 corresponds to equation ( 13 ) with N cu equal to zero. In this 
form of the equation | does not appear explicitly. Therefore it is possible, 
by separating variables, to obtain an analytic solution and to determine the 
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unknown quantity A appearing in equation (13). With the inclusion of angular 
velocity the variables cannot be separated and arriving at an analytical solu- 
tion is not presently possible. 


Procedure for Obtaining Two-Dimensional Solution Curves 


The procedure used to obtain solutions in the present study was to trans- 
form equation (13) into a system of two first-order equations in the following 
manner: 


Let 



( 14 ) 


then 



(15) 


Equations (lA) and (15) were then solved simultaneously by applying the fourth- 
order Runge-Kutta method of numerical integration on a digital computer. (See 
ref . 10 . ) 

The difficulty in applying the method of numerical integration is the 
selection of input data that will yield a solution curve which satisfies all 
the boundary conditions. As part of the data, a point on the free surface and 
the free surface slope at that point must be specified. If the intercept point 
of the free surface with either the tank wall or the T]-axis is known, the ques- 
tion of initial free surface point and slope is answered. However, in general, 
neither intercept is known. It is necessary, therefore, to guess at a value, 
however refined the technique of guessing may be. In addition, a value for the 
constant A must be assigned. The initial value of A was based on an ana- 
lytic expression in reference 6 for the case of zero rotation and a value of tj 
some place along the T]-axis where the slope is known to be zero. 

Thus sufficient information is known to specify the computer input data 
for a given case. However, the solution curve obtained for this set of Input 
data does not, in general, satisfy either the contact angle or the vapor area 
constraint. Therefore, an iterative procedure was incorporated in the computer 
program to converge upon values of A and the rj-axis intercept point with the 
free surface that would yield the solution curve which satisfied the prescribed 
boundary conditions. 

The Iterative procedure used was to compute the solution curve corresponding 
to the first guesses for A and the r \ -intercept denoted as A 0 ,t) 0 . For each 
such computation the contact angle 8 C and the vapor area ratio V bounded by 
the curve were received as computer output data. Then A q and r \ 0 were 
changed slightly and cases were computed corresponding to (A Q +AA, rj 0 j and 
(Aq^o+^o) • By using the information from these cases, the necessary changes 
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in A Q and 


■Ho were approximated by the simultaneous solution of the equations 


d Q _ 9 c( A o +AA o/ T lo) ~ e c( A o^o) ^ + 0 c ( A o^ t 1 o +At 1 o) " ®c ( A o ,t 1 o ) 


Ak r 




dT lo 


dV 


^( A o + ( A o ?T lo) ~ ( A o> ^o ) + ^( A o> T lo +Ar io) " ^( A o ,T lo) 


AA^ 




dT lr 


where 

d 0 c desired change in computed contact angle 

dV desired change in computed vapor area 

A solution curve was then computed corresponding to ( A Q +dA,T] 0 +dT] 0 J and 
the procedure using the new values of (a 0 ,ti 0 ) repeated until a solution curve 
was obtained which satisfied the prescribed boundary conditions* An example of 
the curves utilized in obtaining the solution is seen in figure 3. 



Figure 3*- Representative results of iterative procedure for converging to appropriate A and t\. 
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Determination of an Empirical Formula for A 

In utilizing this technique it was found that solutions to equations (l4) 
and (15) were highly sensitive to variations in A and the r\ -intercept . Usu- 
ally a deviation of a few percent in one of these quantities from the appropri- 
ate value would cause a divergence in the solution, and thereby terminate the 
computation because of demands for numerical values outside the range of com- 
puter capability. Because of the sensitivity of the solution curves to A 

and ti , it was felt that it would be desirable to obtain empirical equations 
o 

for A and t) as functions of the input data, that is, 0 Q , Ng Q , N cu , 

and V. Such data would provide for a more accurate selection and could lead 
to a better physical insight into the free surface shape. 


In order to investigate variations of A with 0 C , N-g Q , N cu , and V, 
each of these parameters was varied independently to determine its effect. 
These variations were then superimposed to obtain a trial function 

A ( e => N Bo> H cu>' r )- 

By setting N-g 0 and N cu equal to zero, equation ( 13 ) may be reduced to 



( 16 ) 


Equation ( 16 ) describes the free surface as a surface of constant curvature; 
a condition satisfied, in two dimensions, by a circle. The variation of A 
with contact angle can be determined by integrating equation (l4) and evaluating 
at the tank wall intercept. Thus, 


(l + V 2 ) 


5/2 


dg = 




(l + T!’ 2 ) 


, 2 !/ 2 


= Ag 


If this equation is evaluated at the tank wall intercept £ - 1, 


A = 


-cot e c 
esc 0 c 


-cos 9 




Computer solutions showed the free surface to be flat for N QU = 0 and 

9 C =. 90°, and hence a surface of zero curvature. For this case equation ( 13 ) 
reduces to 


A " - N Bo% 


where ti is the ti -intercept of the flat surface bounding the appropriate 

° V N Bo V 

vapor area. Since the surface is flat, ti - — and A =■ - where the 

02 2 


vapor area ratio 


, Vapor area 
Total area’ 
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Computer solutions also showed the 
variation of A with N cu to he linear 
when 0 C , Np 0 , and V were held con- 
stant. However, the slope of this lin- 
ear variation was seen to he a function 
of V denoted by F(V). Plots of the 
variation of A with N cu for various 
values of V are shown in figure 4. 
Figure 5 illustrates the nearly exponen- 
tial change in the slopes of th^se lines 
with vapor area ratio* The solid curve 
in figure 5 represents the least squares 
exponential function which hest fits the 
data obtained. 

A linear combination of the fore- 
going results was formed as a trial 
solution for A: 

A(9 c ,N Bo ,N cu ,V) 

N- 

= - + cos e c + F(V) n cu (IT) 


Figure 4.- Variation of A with centrif- 
ugal number for Bond number of 3 and 
contact angle of 90 ° • 


This trial solution for A was found to 
agree with the computer results, in a 



Vapor area ratio, V 


Figure 5.- Slope of A as a function of centrifugal 
number lines for NgQ and 9 C constant. 

A - - + cos e c + F(V) N cu . 
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number of cases to within 0.1 percent, and hence is considered to he a suffi- 
ciently accurate empirical formula for estimating an initial value of A within 
the following region of investigation: 

0 = ®Bo = 6 
0 ^ N cu < 6 
30° = 9 C = 90° 

0.5 = V ^ 3.5 

Bond numbers and centrifugal numbers were chosen in order to cover a large 
range of tanks with liquids exposed to significant capillary forces. The con- 
tact angle was chosen to correspond with those for wetting fluids. However, 
below a contact angle of 30° limitations of the computer technique prohibited 
any further solutions. 


Determination of Bounds of Free Surface ^-Intercept 

As stated previously, the solutions are sensitive to a variation of the 
^-intercept. Thus, t] had to be determined as a function of the boundary con- 
ditions in order to obtain solutions. 

As shown in the preceding section, the liquid free surface is circular 
for N Bo =. N cu =■ 0. Its curvature is given by 


K = -cos 0 C 

For this case the free surface may be represented by an equation of the form: 


! 2 + (r, + d ) 2 


l' 

cos^O 

c 


(18) 


where d is the displacement of the center of curvature from the origin. Equa- 
tion (l8) was solved to yield the variation of the r) -intercept tj from a flat 

surface intercept with contact angle as: 



V 

2 


sec 0 C 


0.5 tan 0 C - 



2 cos^9 c 
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Figure 6.- Examples of rj-intercept variation 
with contact angle. Total tank area, 4. 


An example of the variation of the 
rj -intercept with contact angle is 
seen in figure 6. Since the free 
surface is flat for the case of 90 ° 
contact angle, its ^-intercept is 
given by 



This solution for Hq^c) is plot- 
ted in figure 7 for 0 = 0 C ~ 90 ° • 
The effect of N Bq on tj o ( 0 c ) is 

also shown. It is seen that with 
the introduction of N Bo , the vari- 
ation of tj with contact angle is 
decreased. Figure 8 shows the 
increase in the ^-intercept with 
N C u for constant N Bo and various 
values of 0 C and V. 


Sufficient data were not obtained to determine an accurate empirical for- 
mula for At) q . However, figures 7 and 8 indicate the effects 

of these parameters. For the region investigated the ^-intercepts of the 
free surfaces always lay between that of the flat surface corresponding to 
0 C = 90 °, = N cu = 0; and the circular surface corresponding to 9 C = 0°, 

H-go = N cu = 0. Therefore, the analytical equation for shown in fig- 
ure 7 gives the bounds on the r|-intercept variation. 



Figure 7«- Increase in free surface ^-intercept with 
decreasing contact angle. N cu =• 0. 
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(b) Vapor area ratio of 2. 



(c) Vapor area ratio of 1. 

Figure 8.- Variation of the rj -intercept with 
centrifugal number for Bond number equal 
to 3* 


Possible Application of Two-Dimensional Methods To 
Obtain the Three-Dimensional Free Surface 

It is believed that the three-dimensional free surface to at least a first- 
order approximation can be determined by applying the two-dimensional technique 
simultaneously in two normal directions and defining the three-dimensional sur- 
face by a lattice of such normal two-dimensional solutions. (See fig. 9 *) 

A possible method of approximating the three-dimensional free surface would 
be to obtain an initial two-dimensional curve by the method already presented. 
The second step would be to obtain a number of two-dimensional curves normal to 
and intercepting the initial curve. This second step presents a new problem In 
specifying A in that the computed curve is required to pass through some spe- 
cific point in the space. An investigation was conducted to determine an 
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Figure 9-- Definition of the three -dimensional 
liquid free surface "by a lattice of 
two-dimensional solution curves. 


empirical formula for A under this 
new constraint; however, because of the 
accumulation of considerable computer 
time, this investigation was conducted 
only for the following region: 

0 = N Bo = 2 

0 = N cu = 2 

50° = 9 C = 90° 

% = 1 



The " vapor area” constraint was removed in order not to overspecify the problem. 

The empirical formula for A was obtained by varying one parameter at a 
time and superimposing these results to form a trial formula. The following 
equation is considered sufficiently accurate to approximate A over the region 
of investigation: 


a (0c- i, Bo» h cu) = ‘"Bo - l-22*ou - 


1 + F l( H Bo) 


+ F 2 (”cu) - F l( N Bo) 


cos 0 


where 


F l( N Bo ) “ I “•<= t “(°- 2 5» Bo ) 

F 2 (N cu ) « 0.02168 + 0.2^6l^W cu - 0.01523N CU 2 
RESULTS AND DISCUSSION 


An equation has been derived that governs the three-dimensional equilibrium 
free surface of a contained liquid and which is valid for all magnitudes of 
inertial forces due to translational and rotational accelerations. It has been 
shown that when the inertial forces are sufficiently large to merit disregarding 
the effect of the capillary forces, the equilibrium liquid free surface is 
defined by the intersection of a quadratic equation (in the space coordinates) 
and the region bounded by the tank walls. However, the primary purpose of this 
report is to illustrate qualitatively the effects of various combinations of 
translational and rotational accelerations and capillary forces on the free 
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surface shape in a region in which the inertial forces due to the accelerations 
are of the same order of magnitude as the capillary forces. 

Figure 10 presents examples of the change in the free surface due to 
changes in contact angle for a Bond number of 3, a centrifugal number of 2, 
and a vapor area ratio of 2. The tank is illustrated with a specific height to 
simplify the presentation, but the free surface shape is independent of the 
amount of vapor area added to the bottom of the tank or the amount of liquid 
area added to the top of the tank. It is seen that when the forces due to 
thrust and rotation are sufficiently small with respect to the capillary forces ^ 
the contact angle has a significant effect on the free surface shape. 

Figure 11 illustrates the flattening effect of increasing for con- 

JdO 

stant values of N cu , 9 C , and V. Actually, the effect of N Bo is small in 

the range of interest. It was found that in the absence of rotation, the free 
surface shape was independent of vapor area. This effect is a result of the 
assumption of constant pressure throughout the vapor. 

Figure 12 shows the variation of the free surface with changes in the 
values of N cu * Th e tank is rotating in the plane of the diagram about the cen- 
ter of the tank base. As stated before, the free surface shape is independent 
of any increase in either fluid or vapor area as long as the position of the 
axis of rotation with respect to the fluid is not changed. 

It is seen that with an increase in rotational rate the free surface tends 
to become circular about the axis of rotation. Figure 12 again illustrates the 
pronounced effect of contact angle on the free surface shape since the deflec- 
tion of the surfaces shown away from the circular trend is due to a required 
90 ° contact angle. 

Figure 13 presents the variation of the free surface with changes in vapor 
area ratio. In this figure the axis of rotation is again at the center of the 
tank base. The increases in vapor area shown are completely analogous to 
shifting the axis of rotation vertically away from the fluid. 

For the values of and 0 c corresponding to the cases shown in 

figure 13, the free surface would be flat for N cu - 0. Thus for this combina- 
tion of and 0 C , it is seen that the effect of rotation on the free sur- 

face shape decreases as the distance to the axis of rotation increases although 
the centrifugal force on the free surface increases. This effect is due to the 
fact seen before that the effect of N cu is to force the free surface to become 

circular about the axis of rotation. Therefore the farther the free surface is 
away from the axis of rotation, the greater the radius of curvature and hence 
the flatter the surface which will satisfy the effect of ^ CVi * Equations (8) 

and (9) represent' the extreme case in which the effects of and N cu 

dominate . 
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e c = 90 

e = 6o° 


C 



= 30° 


Tank "base 


Figure 10.- Variations in free surface due 
to changes in contact angle. Bond num- 
ber, 3; centrifugal number, 2; vapor 
area ratio, 2 ; total tank area, 



Vapor 


Figure 11.- Variations in free surface due 
to changes in Bond number. Centrifugal 
number, 0; contact angle, 60°; vapor 
area ratio, total tank area, 4. 



Vapor 


Figure 12.- Variations in free surface due 
to changes in centrifugal number. Bond 
number, 3j contact angle, 90°; vapor 
area ratio, 1; total tank area, 4. 


Figure 13-- Variations in free surface due 
to changes in vapor area ratio, V. 

Bond number, 3j centrifugal number, 6 ; 
contact angle, 90 °; total tank area, 4 . 



CONCLUSIONS 


This report has presented a general equation for the three-dimensional 
equilibrium liquid free surface under the influence of thrust, rotation, and 
liquid-to-solid contact angle. A closed form solution of this equation has been 
obtained which is valid when capillary forces may be neglected. A method has 
been presented for determining solutions to the two-dimensional form of this 
equation by numerical integration on a digital computer. Methods have been pre- 
sented for estimation of the input data required to perform this computer solu- 
tion. Trends in the two-dimensional free surface shape with changes in the 
imposed forces have been shown. 

This investigation was carried out for a specific direction of tank accel- 
eration, position of axis of rotation, and tank geometry. However, the free 
surface equation is valid for the general case. Therefore, various directions 
of thrust, positions of axis of rotation, and tank geometries may be investi- 
gated by appropriate selection of the boundary conditions under which solutions 
to the free surface equation are obtained. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Admini s tr at i on, 
Va., May 25, 1964. 
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